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Preface 


This book is based on a lecture on th^' solution of equations in 
integers which I delivered in 1951 for the participants of a Mathematical 
Olympiad arranged by Moscow State University; I am glad to acknowledge 
the assistance rendered me by my former student. Assistant Professor 
N. M. Korobov, who took notes of my lecture and wrote the first two 
sections and part of the third section pi the book. 

High-school students will readily understand the subject-matter of the 
book. 

A. Gelfond 



Introduction 


The theory of numbers is one of the oldest branches of 
mathematics. It is mainly concerned with the arithmetic properties 
of natural numbers, that is, positive integers. 

A most important problem in what is called analytic theory of 
numbers is the problem of the distribution of prime numbers 
in the sequence of natural numbers. (A prime number is any positive 
integer greater than unity and divisible only by itself and, of 
course, by unity.) It concerns the regularity exhibited by prime 
numbers smaller than some number N for large values of N. 

As long ago as the fourth century B. C. Euclid obtained the 
first results in the solution of this problem. He proved that the 
sequence of prime numbers is infinite. The next result was achieved 
in the second half of the 19th century by the great Russian 
mathematician P. L. Chebyshev. 

Another fundamental problem in, or branch of, number theory 
concerns the representation of integers as sums of integers of some 
specified kind, for example, the possibility of representing odd 
numbers as sums of three prime numbers. This problem (the 
Goldbach conjecture) was solved by the great number theorist, 
the Soviet mathematician I. M. Vinogradov, 

This book is devoted to an interesting branch of number theory, 
the solution of equations in integers. 

The solution in integers of algebraic equations in more than one 
unknown with integral coefficients is a most difficult problem in the 
theory of numbers. The most eminent ancient mathematicians such as 
the Greek mathematician Pythagoras (sixth century B. C.) and the 
Alexandrian mathematician Diophantus (second and third centuries 
A. D.), and also the best mathematicians of more recent times such 
as Fermat (in the seventeenth century), Euler and Lagrange 
(in the eighteenth century) devoted much attention to these 
problems. The efforts of many generations of eminent mathematicians 
notwithstanding, this branch of the theory of numbers lacks 
mathematical methods of any generality, unlike the analytic theory 
of numbers in which many diverse problems can be solved by the 
method of- trigonometric sums, due to Vinogradov. 

As yet, a complete solution of equations in integers is possible 
only for equations of the second degree in two unknowns. Note 
that equations of any degree in one unknown are not really 
interesting: their solution in integers might be carried out by a 
finite number of trials. For equations of degree higher than the 
second in two or more unknowns the problem becomes rather 
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complicated. Even the more simple problem of establishing whether 
the number of integral solutions is finite or infinite presents 
extreme difficulties. 

The theoretical importance of equations with integral coefficients 
is quite great as they are closely linked with many 
problems of number theory. Moreover, these equations are sometimes 
encountered in physics and so they are also important in practice. 
Lastly, the elements of the theory of equations with integral 
coefficients as presented in this book are suitable for broadening 
the mathematical outlook of high-school students and students of 
pedagogical institutes. 

Certain of the main results in the theory of the solution of 
equations in integers have been given here. Proofr of the theorems 
involved are supplied when they are sufficiently simple. 


§ 1. Equations in One Unknown 

Let us consider a linear equation in onh unknown 

n 1 x + a 0 =0 (1) 

with integral coefficients a s and a 0 . The solution of this equation 

_ fJo_ \ 


is an integer only when a 0 is divisible by a x . Thus equation (1) 
is not always solvable in integers. For instance, equation 3x - 
— 27 = 0 possesses an integral solution x = 9, while equation 
5x + 21 = 0 has no such solution. 

The same is true in the case of equations of degree higher 
than the first. For example, quadratic equation x 2 + x - 2 = 0 
has the integral solutions Xj = 1 and x 2 = — 2, whereas equation 
x 2 — 4x + 2 = 0 is not solvable in integers; its roots x 1-2 = 
= 2 + j/2 are irrational numbers. 

The determination of the integral roots of the nth degree equation 

a„x" + a„-ix n ~ l + ... + a x x + a 0 = 0 (n > 1) (2) 


with integral coefficients is not difficult. Indeed, let x = a be an 
integral root of this equation. Then 

a n a " + a„- 1 a"~ 1 + ...+ a x n + a 0 = 0 
a 0 = - a 1 + a„_ 1 a n “ 2 + . . . + 

The latter equality means that a 0 is divisible by a. Consequently 
any integral root of equation (2) is a divisor of its free term a 0 . 
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Check all the divisors of ab one by one; those which transform 
equation (2) into an identity are the integral solutions sought. 
For example, the divisors of the free term of equation 

x 10 + x 7 + 2x 3 + 2 = 0 

are 1, — 1, 2 and — 2. Only one divisor, - 1, is a root of 
the equation; hence this equation possesses only one integral 
root x = - 1. 

Applying the same method it is easy to show that equation 
x 6 - x 5 + 3 x 4 + x 2 - x + 3 = 0 


has no integral solutions. 

The solution in integers of equations in several unknowns 
is much more interesting. 

§ 2. Linear Equations in Two Unknowns 

Let us consider a linear equation in two unknowns 

ax + by + c = 0 (3) 

where a and b are non-zero integers and c is an arbitrary 
integer. We shall suppose that the coefficients a and b have no 
common divisors (except, of course, unity)*. Indeed, if the greatest 
common divisor of these coefficients, d = (a, b), is not unity, then 
a = aid, b = b^, and equation (3) may be written as 


{diX + b x y)d + c = 0 


It can have integral solutions only if c is divisible by d. 
In other words, in the case when (a, b) = d 1, all the coefficients 
of equation (3) must be divisible by d. Cancelling d from the 
equation, we arrive at equation 


diX + bty + Ci = 0 



whose coefficients a y and b x are relatively prime. 

We shall first consider the case c = 0. Equation (3) becomes 

ax + by = 0 (3') 


* Such numbers as a and b are called relatively prime integers. We 
shall denote the greatest common divisor of a and b by (a, b). For relatively 
prime numbers a and b, we have (a, b) = 1. 
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Solving it with respect to x, we obtain 


b 



Obviously, x will be an integer if and only if y is divisible 
by a, or, in other words, if y is a multiple of a, 

y = at 

where t is an arbitrary integer (f = 0, ±1, ±2, . . .). Substituting 
this value of y in the previous equation, we obtain 

b , 

x = at — — bt 

a 

and formulae 

x = — bt, y = at (t = 0, ±1, ±2, ...) 

furnish all the integral solutions of equation (3'). 

We now consider the case c ^ 0. Let us show first of all that in 
order to find all the integral solutions of equation (3), it is 
sufficient to find any one solution, i. e. it is sufficient to find 
integers x 0 , y 0 for which 

ux 0 + by 0 + c - 0 

THEOREM 1. Let a and b be relatively prime and suppose 

[x 0 , y 0 ] ^ tmy solution* of equation 

ax + by + c = 0 (3) 

Then formulae 

x — Xq bt, y = y 0 + at (4) 

where t = 0, ±1, ±2, ..., yield all the solutions of equation 

(3). 

Proof Let [x, y] be an arbitrary solution of equation (3). 

Then equalities 

ax + by + c = 0, ax 0 + by 0 + c = 0 

render 

■ . „ a(x 0 — x) 

ax - ax 0 + by- by 0 = 0; y - y 0 = - y — - 

* A pair of integers x and y which satisfy the equation will be called 
its solution and denoted by [x, y]. 
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Since y - y 0 is an integer and a and b are relatively prime, x 0 — x 
must be divisible by b, i. e. x 0 — * has the form 

x 0 — x = bt 

where t is an integer. But then 

abt 

y - y° = -£-= at 

and we get 

x-x 0 — bt, y = y 0 + at 

Thus, it is proved that each solution [x, y] has the form 
presented as (4). It remains for us to check that any pair of numbers 
[xj, yi] obtained by formulae (4) for an integer t = will be a 
solution of equation (3). To do this, substitute Xi = x 0 — bt u 
y\ — J'o + at \ into the left-hand side of equation (3): 

axi + byi + c = ax 0 — abt x + by 0 + abt l + c = ax 0 + by 0 + c 

Now [x 0 , yo] is a solution, and so ax 0 + by 0 + c = 0 and, 
consequently, 

axi + by t + c = 0 

i. e. [xj, yj] is a solution of equation (3). The proof of the 
theorem is now complete. 

Hence, if one solution of equation ax + by + c — 0 is known, all 
the other solutions can be determined from arithmetic progressions 
whose general terms are 

x = x 0 btj y = y 0 + at (t = 0, ± 1, ± 2, . . .) 

Note that for the case c = 0 the formulae found previously 
x = — bt, y = at 

may be derived from the formulae just derived by setting 
*o = yo = 0- This is legitimate because the values x = 0, y = 0 
are a solution of equation 

ax + by — 0 

But how is one to find a solution [x 0 , y 0 ] of equation (3) 
in the general case when c ^ 0? Consider an equation 

127x- 52 y+ 1 =0 

Let us transform the ratio of the coefficients of the unknowns, 
beginning by isolating the integral part of the improper fraction 


H 



127/52: 


127 „ 23 

ir* 2+ a 


The proper fraction 23/52 is of course equal to 


52/23 


, so that 


127 

52 


= 2 + 


52/23 


Now we shall apply the same transformation to the improper 
fraction 52/23 in the denominator of the last equality: 

52 6 1 

— = 2 H = 2 + — 

23 23 23/6 

The initial fraction is thus equal to 


127 

52 


= 2 + 


2 + 


23/6 


Again, let us repeat the same process for the fraction 23/6: 


23 , 5 . 1 

_- = 3 + _ = 3 + _ 


Then 


127 

52 


= 2 + 


2 +‘ 


3 + 


6/5 


Lastly, we shall isolate the integral part of the improper 
fraction 6/5: 

6 ‘ 1 

T = 1+ T 


The final result is 
127 


52 


= 2 + 


2 + 


3 + 


1 


1 + 1 ~ 
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This expression is a terminating continued fraction. If we omit 
the last term, the one-fifth, and transform the new continued 
fraction so obtained into a common fraction and subtract it 
from the initial fraction 127/52, we get 


2 + 


1 


2 + 

127 

52 


1 


22 

9 


= 2 + 


1 



1143 - 1144 
52-9 



1 

52-9 


22 

9 


Reducing this expression to a common denominator and rejecting 
it we get 


127-9 - 52-22 + 1 =0 


A comparison of this equality with equation 
127 x — 52 y + 1=0 


shows that x = 9, y = 22 is a solution of the equation and by 
the theorem all its solutions are contained in the arithmetic 
progressions 

x = 9 + 52 f, y = 22 + 127 1 (f=0, ±1, ±2,...) 

The result obtained suggests that in the general case of the 
equation ax + by + c = 0 the solution may also be derived by 
expanding the ratio of the coefficients of the unknowns as a continued 
fraction, omitting the last term, and continuing calculations as we 
did above. 

To prove this supposition we shall need some properties of 
continued fractions. Consider an irreducible fraction a/b. Let us 
divide a by b and denote the quotient by q x and the remainder 
by r 2 : 

a = q x b + r 2 , r 2 < b 

Now divide b by r 2 , denoting the quotient by q 2 and the remainder 
by r 3 . Then 

b = q 2 r 2 + r 3 , r 3 < r 2 
Continuing this process we obtain 

r 2 = q 3 r 3 + U, r 4 < r 3 


3-234 
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r 3 = QaU +■ r 5 , r s < r 4 


The quantities q t , q 2 , ... are called partial quotients, and the 
process of calculating them just described is known as the 
Euclidean algorithm. As noted above, the remainders r 2 , r 3 , . . . , satisfy 
inequalities 

b> r 2 > r 3 > r 4 ... ^ 0 (5) 


and thus constitute a sequence of decreasing nonnegative numbers. 

Since the number of nonnegative integers which do not exceed b 
cannot be infinite, the remainder r will vanish at some step and 
the process of forming the partial quotients will cease. Let r„ 
be the last non-zero remainder in sequence (5). Then r n+ , = 0 
and the Euclidean algorithm for the numbers a and b will be 


a = q 3 b + r 2 
b = q 2 r 2 + r 3 

r 2 = q 3 r 3 + U 


( 6 ) 


r «-2 = ‘Zn-i'Vi + r„ 
r„-i= q n r„ 

Let us write these equalities in the form 
a 1 

T- q ' + JK 

b 1 

r 2 ~ q2+ r 2 /r 3 


r n - 2 

r„-i 


— q n - i + 


r n - 1 /r„ 


r rt - 1 


= q n 


By substituting the expression for b/r 2 from the second equation 
into the first equality and the expression for r 2 /r 3 from the 
third equation (which is not written out above) into the second 
equality and so on, we obtain the expansion of a/b as a continued 
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fraction: 


1 


a 


T = qi 


+ 


<h + 


1 

93 + ■ 


+ 


1 


g»-i + 


9 n 


The expression obtained by omitting all terms of a continued 
fraction starting with some particular term is called a convergent. 
The first convergent §! is obtained by omitting all terms 
starting with l/q 2 : 


8i 



the second convergent, 5 2 , is obtained by omitting all terms starting 
with i/q 3 : 


§2 



Similarly, 


83 = 9 i + 


1 


92 + 


93 



84 = 9 i + 


92 + 


93 + 


94 



etc. 

Because of their method of formation the convergents satisfy 
the obvious inequalities 

81 <83 < . . . < §2i- 1 < -g 

8 2 > 64 > . . . > 5 2 * > ~ 
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Let us write the kth convergent 8* as a fraction 
5* = ^ (1 ^ k < n) 


and find the rule for forming the numerators and denominators 
of convergents. We begin with the first three conver gents 8 1( 8 2 
and S 3 : 

81 = 9i = = ~q^’ = 9i> Qi = 1 

6 2 = 9l + P 2 = qi q 2 + 1; G 2 =q 2 
92 92 «2 


83 — 9i + 


92 + 


1_ 

93 


- 9i + 


93 


9293 + 1 


9i9293 + 9i + 93 
9293 + 1 


P 3 = 9i 929 3 + 9i + 9s; 63 = 9293 + 1 
From these we obtain 


P3. 

63 


P3 — P 293 + ^*1 » 63 ~ 6293 + 6l 

By applying mathematical induction * we can prove that the 
similar relations 


Pk — Pk-lQk + ^*-2. Qk ~ 6*-l9* + 6*-2 (7) 

hold for all k ^ 3. 

Let equalities (7) be valid for some k ^ 3. By the definition of 
the convergents it immediately follows that if in the expression for 

8*, q k is replaced by q k + — - — , 8* becomes 8*+ 1 . By the induction 

9 *+ 1 


hypothesis. 


P k _ Pk- l9* + Pk-2 
Qk Qk-iVk + Qk-2 


The substitution of q k by q k H 

9 *+ 1 

changes the latter to 8 k+1 so that 


in the expression for 


5 , 


* See I. S. Sominsky, The Method of Mathematical Induction, Mir 
Publishers, Moscow. 
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§*+l — 


p k + —~p k -1 
4*+ 1 

6* + — — Qk-i 
4*+l 

•P*4*+l + Pk-1 

6*4*+ 1 + Qk-i 


Hence, as 5* +1 = 


it follows that 


^*+1 — ^*4k + ?k~U 6*+i - 6*4*+ 1 + 6*-l 

Thus, if equalities (7) hold for some k > 3 they are also valid 
for k + 1. For k = 3 equalities (7) are indeed satisfied, so they are 
valid for every k ^ 3. 

Let us now show that the difference between consecutive convergents 
8 k — 8*_ j satisfies the relation 

8 *- 5 -'*S £ r (t>,) (8) 


Indeed 


<• _ s _ P k P k- 1 _ PkQk - 1 ~ QkPk- 1 

k 6* G*-1 6*6*- 1 


Using formulae (7) we can transform the numerator of this fraction 

PkQk- 1 - QkPk- 1 = (Pk- l4* + Pk-i)Qk- 1 — 

— (6*-l4* + 6*— 2 ) Pk— 1 = ~ (Pk-lQk -2 — Qk-iP *- 2 ) 

The expression in brackets is obtained from the initial one by 
replacing k by k — 1. 

Repeating similar transformations we get a chain of equalities 

PkQk- 1 ~ QkPk-1 l)(^ > *-l6*-2 ~ 6*- \Pk-l) = 

= ( — l) 2 (Pk-2Qk - 3 — Qk-2Pk-i) — ■■■ 

••• = (— l) t_2 (F 261 — Q 2 P l) = 

= (- l) k ~ 2 (4i42 + 1 - 424 i) = (- l)* -2 
whence it follows that 

s _ s PkQk-k-QkPk-k ^ (-I)*' 2 (-1 f 

k k 1 6*6*- 1 6*6*- 1 6*6*- 1 





If the expansion of a/b into a continued fraction contains n 
terms, then the nth convergent 8, will coincide with a/b. Applying 
equality (8) for k = n, we get 


- 8,-1 = 


(- 1 )" 

6b6b-1 



(-1)" 
bQ n - 1 


(9) 


We return now to the solution of equation 

ax + by + c = 0, (a, b)= 1 (10) 


We rewrite relation (9) in the form 

fl Pn - 1 (-1 r 

8 6,- 1 bQ„-i 


Reducing the fractions to a common denominator and discarding 
it we obtain 


^6a— i bP n ~ i — ( 1)° 

a6n-l+b(-J > a-l) + (-ir 1 =0 


Let us multiply this expression by (— 1)" 1 c. Then 

<* [(— 1)“ ‘ c 6n- i] + 8 [( — 1)" cP n ~ i] + c = 0 
Hence the pair of numbers [x 0 , y 0 ]. suc ^ that 

x 0 = (-ir 1 c6„-i, yo = (— lfcP,-! (li) 

is a solution of equation (10); according to theorem 1, all solutions 
of this equation are of the form 

x = {-lY~ 1 cQ n . 1 -bt, y = (- lfcP,-, +at 

(r = 0, ±1, ± 2 , ...) 

This fully solves the problem of determining all the integral 
solutions of linear equations in two unknowns. 


§ 3. Equations of the Second Degree in Three 
Unknowns (Examples) 


EXAMPLE 1. Consider equation 

x 2 + y 1 — z 2 (12) 

From a geometrical point of view, the determination of integral 
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solutions of this equation amounts to finding all Pythagorean 
triangles, i. e. right triangles whose legs x, y and hypotenuse z are 
represented by integers. 

Let us denote the greatest common divisor of the numbers x and y 
by d: d = (x, y). Then 

x = Xid, y = y x d 
and equation (12) becomes 

x\ d 2 + y\d 2 = z 2 

This means that z 2 is divisible by d 2 , and hence z is a multiple 
of d: z = Zjd. 

Equation (12) can now be written as 

x 2 d 2 + yjd 2 = z 2 d 2 

Cancelling d 2 we get 

xi + yl = z\ 

This is an equation of the same type as the initial one, (12), 
only x t and y x have no common divisors (except, of course, 
unity). So when solving equation (12) we can restrict ourselves 
to the case when x and y are relatively prime. 

Thus we may suppose that (x, y) = 1. Then at least one of the 
quantities, x and y (say x), is odd. Transferring y 2 into the right-hand 
side of equation (12) we get 

X x 2 = z 2 — y 2 , x 2 = {z + y){z - y) (13) 

We shall denote the greatest common divisor of the expressions 
z + y and z — y by d v Then 

z + y ad u z - y -bdi (14) 

where a and b are relatively prime. Now substituting the values 
of z + y and z — y into (13) we obtain 

x 2 = abd 2 

Since a and b possess no common divisors, the latter equality 
is possible only if these numbers are perfect squares*: 

a = u 2 , b = v 2 


* The product of two relatively prime numbers is a perfect 
square only if each factor is a perfect square. 
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But then 

x 2 = u 2 v 2 d\ 

and 

x = uvd x (15) 

Now determine y and z from equalities (14). Adding them together 
we get 

2z = ad x + bdi — u 2 d x 4- v 2 d 1 ; z — U + V d x (16) 

while subtracting the second one of equations (14) from the 
first we get 

u 2 - v 2 

2y = ad x — bd t = u 2 d x — v 2 d x ; y = d x . (17) 

From (15) it follows that, x being odd, u, v and d x are also odd. 
Moreover, d x — 1, since otherwise from the equations 

j A u 2 v 2 j 
x = uvd x and y = — — d x 

it would follow that x and y have a common divisor d x # 1, 
which contradicts the supposition that they are relatively prime. 
The numbers u and v are connected with the relatively prime 
numbers a and b by the equations 

a 2 = u 2 , b = v 2 


and so are relatively prime themselves; v < u since b < a, as can be 
seen from (14). 

Substituting d x = 1 into equalities (15)-(17) we get formulae 


x = uv, y = 


tr + v 2 


which, with odd and relatively prime u and v (v < u), furnish all the 
triplets of positive integers x, y, z which do not possess common 
divisors and which satisfy equation (12). By a substitution of the 
expressions for x, y and z in equation (12), it is easy to verify 
that for arbitrary u and v the numbers (18) satisfy this equation. 

For the initial values of relatively prime u and v, formulae (18) 
yield the following frequently encountered equalities 

3 2 + 4 2 = 5 2 (v = 1, u = 3) 

5 2 + 12 2 = 13 2 (v = 1, u = 5) 
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15 2 + 8 2 = 17 2 (t> = 3, u = 5) 

As was noted above, formulae (18) give only those solutions of 
equation 

x 2 + y 2 = z 2 

in which the numbers x, y and z do not have common divisors. 
All the rest positive integral solutions of this equation can be 
obtained by multiplying solutions (18) by an arbitrary common 
factor d. The method used for determining all the solutions of 
equation (12) can also be employed to find all the solutions of 
other equations of the same type. 

EXAMPLE 2. Find all the positive integral solutions of equation 

x 2 + 2 y 2 = z 2 (19) 

if the numbers x, y and z are pairwise relatively prime. 

Note that if the triplet x, y, z is a solution of equation (19) 
and the numbers x, y and z possess no common divisors 
(except, of course, unity), then they are pairwise relatively prime. 
Indeed, let x and y be multiples of a prime number p (p > 2). 
Then from equality 



with an integral left-hand side it follows that z is a multiple of p. 
The same conclusion holds if x and z, or y and z are multiples of p. 

Notice that x must be an odd number for the greatest common 
divisor of x, y and z to be equal to unity. For if x is even, then the 
left-hand side of equation (19) is an even number so that z is 
also even. But then x 2 and z 2 are multiples of 4. From this it 
follows that 2y 2 is divisible by 4, in other words that y must also 
be an even number. Thus, if x is even then all three numbers 
x, y , z must be even. Thus, in a solution not having a common 
divisor different from unity x must be odd. From this it immediately 
follows that z must also be odd. Transferring x 2 into the right-hand 
side of equation (19) we get 

1y 2 = z 2 - x 2 = (z + x)(z - x) 

But z + x and z — x have the greatest common divisor 2. Let 
their greatest common divisor be d. Then 

z + x = kd, z — x = Id 

where k and / are integers. Adding together these equalities, 
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then subtracting the second one from the first we arrive at 
2 z = d (k + /), 2 x = d (k — l) 

But z and x are odd and relatively prime. Therefore the greatest 
common divisor of 2x and 2z must be equal to 2, that is d = 2. 

z + x z — X 

Thus, either — - — or — — — is odd. Therefore either 
2 2 

, z - x 
z + x and — - — 


are relatively prime or 


z + x 


and z — x 


are relatively prime. 

In the first case equality 

(2 + x) —j— = y 

leads to 

z + x — n 2 , z — x = 2 m 2 
while in the second case from 


z + x 


(z — x) = y 2 


it follows that 


z + x = 2 m 2 , z — x = n 2 


where n and m are positive integers and m is odd. Solving these two 
systems of equations with respect to x and z, and finding y, 
we obtain either 

z = -y(” 2 + 2m 2 ), x = y(n 2 -2m 2 ), y = mn 
or 

z = y(n 2 + 2m 2 ), x = -^-(2 m 2 - n 2 ), y = mn 

respectively, where m is odd. Combining these two expressions 
we derive the general formulae 


1 


1 


x = ± -y(n 2 - 2m 2 ), y = mn, z = —(n 2 + 2m 2 ) 
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where m is odd. But for z and x to be integers, n must be even. 
Putting n = 2b and m = a, we finally obtain general formulae which 
yield all the solutions of equation (19) in positive integers x, y and z 
having no common divisors greater than unity: 

x = ±(a 2 — 2b 2 ), y — lab, z = a 2 + 2 b 2 (19') 

where a and b are positive and relatively prime and a is odd. 
No other restrictions are imposed on a and b except that x 
should be positive. Formulae (19') do indeed provide all the 
solutions in integral and relatively prime x, y and z, since on the one 
hand we have proved that in this case x, y, z must be 
represented by formulae (19'), while, on the other hand, any 
numbers a and b complying with the conditions formulated above 
furnish such relatively prime numbers x, y, z as constitute a 
solution of equation (19). 

§ 4. Equations of the Type x 2 — Ay 2 — 1. 

Finding All Solutions of This Equation 

We now come to the solution in integers of equations of 
the second degree in two unknowns of the type 

x 2 - Ay 2 = 1 (20) 

where A is a positive integer other than a perfect square. To find an 
approach to the solution of such equations, let us expand 
irrational numbers such as \/a into continued fractions. From 
Euclid’s algorithm it follows that any rational number may be 
expanded into a continued fraction with a finite number of terms. 
For irrational numbers the situation is different: their expansions 
into continued fractions are infinite. 

Let Us find, for example, the continued fraction expansion 

of the irrational number |/2. Consider an obvious identity 

0/2 - l)(]/2+ 1) = 1 
or 


]/2+l 


2 + (|/2 - 1 ) 

Replacing the difference ]/2 — 1 in the denominator of the last 
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identity by the expression 


2 + 0 - 1 ) 


which is obviously equal to it, we receive 


]/ 2 - 1 = 


l/2 = i + 


2 -+ ( 1/2 — 1 ) 


2 + (j/2 — 1) 


Again we replace the bracketed term, in the denominator of the last 
equation, by the fraction equal to it from the same identity. Then 
. . 1 


l/2=l+- 


2 + (l/2 — 1) 

Continuing this process, we arrive at the following expansion of 
j/2 into an infinite continued fraction 

j/2=l+ - (21) 


Note that the method of expansion based on identities of the type 
(1 /m 2 + 1 - m)(]/m 2 + 1 + m) = 1 

is not suitable for all irrational numbers ]/a. It may obviously 
be used when the integer A may be expressed as A = m 1 + 1 
where m is a non-zero integer. (In particular, the case m— 1 
leads to the expansion for + = ]/l, m = 2 corresponds to A = j/5 etc.) 
However, comparatively simple methods also exist for the expansion 
of j/a into continued fractions in the general case. 

As before, in the case of finite continued fractions, we shall 
form for the infinite continued fraction (21) a sequence of convergents 

Si, 82, 83, ... 

81 = 1, 8j <l/2 
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5 2 = l+y = 4’ 52> J^ 

5 3 = 1 + — - — = - j , 63 < ]/2 (22) 

2 + l 

84 = • = iy’ 84 


etc. 

From the way these convergents are formed it follows that 

8 t < 63 < ... < j/2 
6 2 > 84. > . . . > |/2 

In general, if we are given the continued fraction expansion of some 
irrational number a 


a = 4, + 


1 


<?2 + 


1 

«3 + 


then the convergents satisfy the inequalities 

81 < 83 < . . . < 8 2ik+ 1 < . . . < a < . . . ^23) 

. . . < & 2k < ■ ■ ■ < 84 < 8 2 
Let us write the convergent 8* as 

s Pk 

‘-Q* 

Expressions (7) 

Pk — Pk-lQk + P *- 2 > Qk = Qk-lVk + Qk - 2 

derived in § 2 for the case of finite continued fractions are also 
valid for infinite fractions, as in the derivation of (7) we did not 
make use of the fact that the continued fraction was finite. 
Hence relation (8) between consecutive convergents 

* * _ ("IF 


■ i — 


QkQk- 1 


also remains valid. 
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5 3 — §2 — 


Assume for example fej = 3 and k 2 = 4 and expand ]/2 into a 
continued fraction. Equalities (22) will then lead to 

7 3 -1 

5 2 ~ 10 

c s _ 17 7 1 

5 4 -6 3 - — - T = — 

which coincides with the results given by formula (24). 

Consider now formula (24) for the subscript 2k: 

(- l) 2l+1 1 


§2* — 8 2t+1 - — (5 2)[+ i — 5 2J[ ) — — n n 

y2*+lV2k Y.2k+l\l2k 

We shall now prove the validity of inequality 


0 < P 2k ~ U-Q.2k < 


1 


(25) 


62*+ 1 

The left inequality is obvious, for, according to inequalities (23), 

p 

a < 8 2 * = ; <*62* < E 2jt ; 0 < P 2k — v-Qik 

The deduction of the other inequality (25) is also a rather simple 
procedure. From (23), 

§2*+ 1 < a < S 2lt 


so 


5 2 * — < b 2k — 8 2k + x — 


1 


62*62*4 


Substituting P 2 k/Q 2 k for 8 2k , we get 


2* 


- a < 


1 


62* 62*62*4-1 

Multiplying this inequality by Q 2k we arrive at the desired 
result: 


P 2 k — <*62* < 


1 


62*4- 1 

We now apply the results obtained to the solution of equation 

x 2 — 2y 2 = 1 (26) 

Let us transform the left-hand side of this equation: 
x 2 - 2y 2 = (x - |/2y)(x + \/2y) 
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and assume * = P 2k and y = Q 2k , where P 2k and Q 2k are the 
numerator and the denominator, respectively, of the corresponding 
convergent in the expansion of ]/l. Then 

P\ k - 221* = (P 2k - l/2 Q 2k ) (p 2k + 1/2 Q 2k ) (27) 

The left-hand side of this equality, and therefore the right- 
hand side too, is an integer. We shall show that this integer is 
greater than zero but less than two and so is equal to unity. 
To do this write inequality (25) for a = (/2: 

Q < P 2k — l/2 Q 2k < -yr— — • . (28) 

V. 2*+l 

From this it is clear that both factors of the right-hand side of 
(27) are positive, and so 

P 2k — 2 2l* > 0 


On the other hand. 


P 2k 1/2 Q 2k < 


1 

Q 2k +i 


1 


1 


l 


Q 2k y 2k +i + Q 2k -i 2 Q 2k -i- Q 2k — 1 222 * 


But, because of inequalities (23), 




Hence 

l/222* < P 2k 
P 2k |/2 Q 2k < 2 P 2k 

and the factors on the right-hand side of (27) satisfy inequalities 

P 2k l/222* < 

^<Z2k 


P 2k l/2 22* < 2 P 2k 

Multiplying these inequalities together gives 
^21 ~ 22l* < 


Q 2k 


Using inequality (28) we arrive at 
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P\k ~ 2 Qi* < 
For any k > 1 


-g i 

Qlk Q.2kQ.2k+l 


1 1 

< 


therefore 


Q.2kQ-2k + 1 6283 10 

Plk-2Qlk<]/2+-^<2 


We have thus proved that for any k ^ 1 the integer Pj k — 
— 2 Ql k satisfies inequalities 

0 < P\k — 2 Qlk < 2 

Hence 

Pik-2Ql k = l 

This means that for any k > i the numbers x — P2h y — 0,2k 
yield the solution of the equation 

x 2 — 2y 2 = 1 

We do not yet know whether or not the solutions of equation 
(26) found above are all the solutions of that equation. 

The question now naturally arises, how do we find all the solutions 
of equation 

x 2 -Ay 2 = l (29) 

in integers x and y for integral A > 0 and irrational j/A? 
We shall show that we can do this if we can find at least one 
solution of equation (29). As evidenced by equation (26) such 
equations do have solutions. So we shall now consider the problem 
of how to obtain all the solutions of equation (29) from a single 
particular solution which we shall call a minimum or least solution 
leaving open for the moment the question of whether or not 
equation (29) always has at least one solution in integers other 
than the trivial solution x = 1, y = 0. 

Let us suppose that equation (29) does have a non-trivial solution 
[*o, J’o]. *0 > 0, y 0 > 0, and 

x o - Ayl = 1 (30) 

(Remember that a solution is a pair of integers [x 0 , y 0 ] satisfying 
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the equation.) We shall call this solution minimal if for x = x 0 
and y=y 0 the binomial x + J [a y, ][~A > 0, assumes the least 
possible value among all the possible values which it will take 
when all the possible positive integral solutions of equation (29) 
are substituted for x and y. For example, the least solution of 
equation (26) is x = 3, y = 2 because for these values of x and y 

the binomial x + ]/ly assumes the value 3 + 2 ]/l. Indeed, equa- 
tion (26) admits of no other solutions with small positive integers 
x and y; the smallest values of x and y constitute the next 
solution: x = 17, +=12 and it is clear that 17 + 12J/2 is 
greater than 3 + 2 j/2. Note that equation (29) does not have two 
least solutions. For, assume that solutions [xi, y j] and [x 2 , y 2 ] 
give the same value to the binomial x + y A y. Then 

*i +Y^yi =x 2 + ]/Ay 2 (31) 

However, \f~A is an irrational number while x lf y t , x 2 , y 2 are 
integers. Hence, as it immediately follows from Eq. (31) 

~x 2 =(y 2 ~yi)\/A 

which is impossible because x x — x 2 is an integer and (y 2 — +) ]/~A, 
being a product of an integer and an irrational number, is irrational. 
And we know that an integer cannot be irrational. The contradic- 
tion disappears if x x = x 2 and y 2 = y 2 , i. e. if we take not two different 
solutions, but one. Thus, if a least solution does exist, it is unique. 

Observe now another very important property of the solutions of 
equation (29). Let [xj, +] be a solution of this equation. Then 

x\ - Ay\ = 1 
or 

(*i + l/iyi)(xi - j/iyO = 1 (32) 

Now raise both terms of equality (32) to the positive integral 
power n: 

(x 1 +l/5yi)”(x 1 -lA4y 1 ) n = l (33) 

Raising the factor on the left-hand side to the power n according 
to the binomial theorem, we get 

(xi + l^AyiT = A + nxT 1 ]/Ay 2 + 

+ n(n-0_x-- 2i 4 y f + ... + ( ]/Aryl = x m + ]fAy n (34) 

2 

where x„ and y n will be integers since the first, the third term, and, 
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in general, the odd terms of the binomial expansion are integers 
while the even terms are integers multiplied by j/A. Collecting 
separately the odd and the even terms of the expansion we 
obtain (34). We shall now prove that the numbers x„ and y„ 
will also be a solution of equation (29). The proof is simple: changing 
the sign of ]/~A in equality (34), we obtain 

(*1 -• V^yiY = *. - ]/Ay„ (35) 

Multiplying (34) and (35) term by term and using expression 
(33) we finally have 

(*i + V^yiYixi - V^yiT = 

= (x„ + |/Ay„)(x„ - ]/Ay n ) = x\ - Ay\ = 1 (36) 

or, in other words, [x„, y„] is also a solution of equation (29). 

Now we can prove the basic theorem concerning solutions 
of equation (29): 

THEOREM II. Any solution of equation (29) 
x 2 - Ay 2 = 1 

with positive A and irrational ]/a is of the form [± x n , ± y„] where 

X„=~ [(*0 + >’o V A f + ( x 0 - y<) 1/A)"] 

y„ = ~=- [(x 0 + yoV A ) n - (x 0 - yol/A)"] 

2]/A 

and [x 0 , y 0 ] is the least solution of the equation. 

Proof. Suppose the converse, namely, that there exists a positive 
integral solution [x', /] of equation (29) such that the equality 

x' + 1 /Ay' = (x 0 + 1 fAyf)" (38) 

does not hold for any positive integer n. Consider a sequence of 
numbers 

X 0 + ]/~A yo, (x 0 + 1 /a yo) 2 , (Xo + i/a y 0 ) 3 . • • • 

It is a sequence of positive and indefinitely increasing numbers,* 
since x 0 > 1, y 0 > I and x 0 + (/Ay 0 > 1. 

By definition of [x 0 , y 0 ] as the least solution, 

x' + 1 /Ay' > x 0 + l rky 0 

there always exists an integer n > 1 such that 

(x 0 + 1 /Ay of < x' + I /a/ < (x 0 + 1 /Ay 0 )" + 1 (39) 
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But x 0 - ]/Ay 0 > 0 because 

(x 0 + ]fAy 0 )(x 0 - ]/~Ay^) = xl - Ayl = 1 > 0 

Hence, when all terms of inequalities (39) are multiplied by the 
same positive number (x 0 - \/a y 0 ) n the inequality signs are 
retained, and we will have 

(x 0 + 1 /Ay 0 )"(x 0 ~ ]fAyoT < (*' + l/~Ay')(x 0 - ]/Ay 0 ) n < 

< (x 0 + \/Ay 0 ) n+ 1 (x 0 - ]/Ay o y (40) 

Since 

(x 0 + |//ty 0 )" (x 0 - | [Ay 0 ) n = (xg - Ay 2 0 ) n = 1 (41) 

we have 

(x 0 + \/Ay 0 ) n+ 1 (x 0 - l /Ay 0 ) n = x 0 + j/Ay 0 (42) 

In addition to this 

(x' + \ /Ay')(x 0 - ]fAy 0 ) n = (x' + ) /Ay')(x„ - ]/Ay„) = 

= x'x„ - Ay'y„ + ]/a (y'x„ - x'y„) = x + l/ Ay (43) 

where x and y are integers and 

x n - V A y* = ( x o - V^yoT 

Making use of relations (41)-(43) and inequalities (40), we obtain 
inequalities 

1 < x + ]/~Ay < x 0 + 1 /Ay 0 (44) 

We shall show that the pair of integers x and y is a solution of 
equation (29). To do this, multiply termwise, Eq. (43), i.e. equation 

x + \/Ay = (x' + 1 /Ay')(x 0 - j/iyo)” (45) 

and equation 

x - ]/Ay = (x' - ]/ r Ay')(x C) + \/Ay 0 )" (46) 

which is immediately obtained from (43) by a change of the sign 
of ]/a. As [x', /] and [x 0 , y 0 ] are solutions of equation (29), 
the result will be 

(x + | /Ay)(x - l/Ay) = x 2 - Ay 2 = 

= (x' + ]/A /) (x' - 1 /A y') (x 0 + l/A y 0 )" ( x 0 ~ ]/a y 0 T = 

= (x' 2 ^Ay' 2 )(x 2 -Ay 2 0 r=l (47) 
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The last step is to prove that both x and y are positive. First 
of all, note that x # 0, otherwise (47) would give us 

-Ayl = l 

which is impossible because A > 0. Moreover, if y - 0 the same 
equality (47) furnishes x 2 = 1, but inequalities (44) yield x > 1, a 
contradiction. Finally, note that the signs of x and y coincide. 
For if we suppose that the signs of x and y are opposite, 
then those of x and — y are the same. Let us compare the moduli 
of the numbers x + \/~Ay and x — ]/Ay. The modulus of the first 
number must be less than of the second because, in the first 
case, two numbers with the same sign are subtracted one from the 
other while in the second case they are added together. But we 
already know that 

x + ]/~Ay > 1 

and so x — \/Ay is also greater than unity in modulus. But 

(x + ]/Ay)(* — j//4y) = x 2 - Ay 2 = 1 

and we have arrived at a contradiction, because the product of two 
numbers, each with modulus greater than unity, must also be 
greater than unity in absolute value. Hence the signs of the two 
numbers, x and y, are the same and x # 0 and y ^ 0. Inequalities 
(44) lead then directly to the conclusion that x > 0 and y > 0. 
And so, by supposing that there exists a solution [x', y'] of equation 
(29) 

x 2 — Ay 2 = 1, A > 0 

such that (38) does not hold for any positive integer n, we have 
been able to construct another positive integral solution [x, y] of 
the same equation (x > 0, y > 0 are integers) satisfying inequalities 
(44), which contradict the definition of the least solution [x 0 , y 0 ]. 
We have thus proved that the supposition that there exist solutions 
not given by formula (38) leads us to a contradiction. In 
other words we have proved that all solutions of our equation may 
be obtained from formula (38). 

Thus each solution [x, y] of equation (29) may be derived from the 
expression 

x + ]/Ay = (x 0 + ]/Ay 0 f , n^0 (48) 

where [x 0 , y 0 ] is the least solution. Changing the sign of 
J [a in equality (48), we also get equality 

x - | [Ay = (x 0 - l/W (49) 
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Adding the two equalities, then subtracting the latter from the former, 
and finally dividing the sum by 2 and the difference by 
2 ]/A, we obtain 


* = = 4- [(*<> + V^yoY + (*o - V^yof] 


y = y„ = 


2 ]/a 


[(x 0 + J /Ay 0 Y - (*o - V^yoYl 


(50) 


These are explicit expressions for any positive solution [x, y]. 
Each solution can be derived from them by an arbitrary choice 
of the signs of x„ and y„. 

For example, since we saw above that the least solution of 
equation x 2 — 2y 2 = 1 is x = 3, y = 2, all the solutions of this 
equation are contained in formulae 


x„ = i-[(3 + 2l/2 Y+ 0-2l/2r] 


* = j^[V + 2\/2Y~(3-2\/2n 

The least solution [3, 2] corresponds to n = 1 ; for n = 2 and 3 
the solutions are [17, 12] and [99, 70] respectively, etc. 

Note that the numbers x„ and y n increase with n as terms of a 
geometrical progression with a common ratio x 0 + y Ay a . Indeed 

0 < x 0 - ]/Ay 0 < 1 

because 

(x 0 + ]/Ay 0 )(x 0 -]/Ay 0 )= 1. 

Therefore, (x 0 — ]/a y 0 )" tends to zero with increasing n. 

Note also that if equation (29) possesses at least one non- 
trivial solution, or, in other words, at least one solution with 

y 0, then it will also have a least solution and, consequently, 

all its solutions can be derived from formulae (50). The existence 

of a non-trivial solution for any positive integer A with 
J /~A irrational will be discussed in § 5. 

§ 5. Equations of the Second Degree in Two Unknowns: 
the General Case 

We shall prove, in this section that for an arbitrary positive 
integer A with ]/A irrational, the equation 

x 2 — Ay 2 = 1 


(51) 
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always has a non-trivial solution; in other words, there exists a pair 
of non-zero integers x 0 and y 0 satisfying the equation. 
We shall first describe a method for expanding an arbitrary 
positive number into a continued fraction. (Previously, when 
expanding ]/l into a. continued fraction (see § 4), we made use 
of the special properties of that number.) Let a be any 
positive number. Then there always exists an integer less than or 
equal to a, and greater than a - 1. This integer is called the 
integral part of a and is denoted by [a]. The difference 
between a and its integral part is called the fractional part 
of a and is designated by {a}. The relation between these two 
quantities 

«-[«]- {“} 

or 

a = [a] + {a} (52) 

is a direct consequence of their definitions. 

Note also that the fractional part of a number, being the 
difference between a positive number and the greater integer not 
exceeding it, is nonnegative and always less than unity. For 
example, the integral part of 27/5 is 5 and its fractional part is 
2/5; for J/2 the respective numbers are 1 and \/2 — 1 , for 

j/^52, 3 and j/52 - 3, etc. 

The notions just introduced may be used for the expansion of 
positive numbers into continued fractions. Suppose 

[«]=«.. M = ~ 

a l 


Then 




(53) 


As {a} is always less than unity, a t is always greater than unity. 
If a were an integer then its fractional part would be equal 
to zero so that 0 ^ would be infinite and we would have 
a = q v However, as we are discussing the continued fraction 
expansion of irrational numbers, we can leave aside this particular 
case and say that oti is a positive number greater than unity. With 
this number we can proceed as we did with a, and write 

a i = q 2 + ~~t -^- = {<* 1 } 
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Repeating this process, we obtain a series 



4i = l>] 

1 

a l — 42 + > 

a 2 

42 = [ai] 

1 

a 2 = 43 + , 

<*3 

43 = [“ 2 ] 


3 

II * 

r““i 

1 


( 54 ) 


For rational numbers a (or, which is the same, for a = a/b, 
where a and b are positive integers) this sequential calculation 
furnishes the same result as the Euclidean algorithm with integers 
q u q 2 , # 3 , q„, ... being the partial quotients (see formula 

(6) in § 2). Here, as also in § 2, the process must break off. 
On the other hand, when a is irrational this process must be 
infinite. For if a„ were an integer for some n, then a„_j 
would be rational, and so also <x„_ 2 , <x„_ 3 , ... and, lastly, a, 
would be rational. Consecutive substitutions eliminating <x lt a 2 , ... 
..., a,,-! from formulae (54) lead to a continued fraction 

1 

a = 4 i + — 


+ , (55) 

or, since n may be taken arbitrarily large, we can write it 
in the form of an infinite continuous fraction 


a = 4i + 


1 


42 + 


1 

43 + 


' + 


1 

4n + 
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As mentioned in § 4, relation (8) between the convergents 
is not dependent on the finiteness or infinity of the continued 
fraction and so it holds also in this case. From expression (8), 
as we have seen, follows inequality (25) for even convergents. 
This inequality will again be used to prove the existence of a solution 
of equation (51), but the reasoning will be more complicated than 
in the particular case when A — 2. 

THEOREM III. For any positive integer A and irrational J f~A 
equation (51) 

x 2 — Ay 2 = 1 

possesses a non-trivial positive solution [x 0 , y 0 ]- 

Proof. Because of certain complications in the proof of the 
existence of solutions to equation (51) we shall break up 
the proof into several steps. The first step will prove the 
existence of a positive integer k such that equation 

x 2 — Ay 2 = k (56) 

has an infinite number of positive integral solutions. Let us 
consider the binomial x 2 - Ay 2 . We shall replace x and y respectively 
by the numerators and denominators of the consecutive- even 
convergents of the irrational number a = |/A. Then 

z 2„ = Pin ~ AQL = (Pin ~ *Q2n)(P 2 n + «Q 2 n) (57) 

But since 

0 < P 2 n - a Qm < — 

Sl2n + 1 

it directly follows that 

0<P 2n + ctQ 2n = 2a Q 2n + P 2n - a Q 2n < 2 a Q 2n + — - — 

\C2n+ 1 

Let us use the last two inequalities to estimate z 2n . By 
substituting greater quantities for both factors on the right-hand 
side of (57) we get for z 2n the inequality 

0 < z 2 n < ( 2aQ 2n + — ) < 2a + 1 (58)" 

VJ2n+l \ SZ2n+ 1 / 

since Q 2 „ is less than Q 2n+l . If we replace x and y by P 2n 
and Q 2n respectively in the binomial 

z = x 2 — Ay 2 

z will assume an integral positive value. Thus, all the numbers 
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z 2 , z 4 , . . . , z 2n , . . . will be positive integers, none of which 
exceed the same number 2a + 1. But since a = ]/a is irrational, 
its continued fraction is infinite and so the sequence of pairs 
of numbers P 2 „ and Q ln is also infinite. Now since there are not 
more than [2a + 1] integers between 1 and the number 2a + 1 
(which is definite and does not depend on n), the infinite 
sequence of positive integers z 2 , z 4 , . . . , z 2 „, ... is made up of 
a finite number of different terms. In other words, the infinite 
number series z 2 , z 4 , . . . , z 2 „, ... is just the sequence of 
integers 1, 2, 3, ... , [2a + 1] repeated in some way or other and it is 
not even necessary for all these integers to occur in the series. 
Note also that since the quantity of different terms of the 
infinite series z 2 , z 4 , . . . , z 2 „, ... is finite, at least one term 
(one number), k (1 < k < [2a + 1]), is repeated an infinite number 
of times. In other words, among the pairs of numbers 
[^ 2 , Qz], IP*, Q*], ■■■, [Pin, Qln ],-•• there is an infinite 
set of pairs for which z = x 2 — Ay 2 assumes the same value k upon 
substitution of these numbers in place of x and y. Thus, we have 
proved the existence of a positive integer k for which equation 
(56) possesses an infinite number Of integral solutions [x, y\ Let 
us enumerate once again these number pairs which are solutions 
of equation (56) for given k denoting them by [u 1; t^], [u 2 , t> 2 ], . . . 

. . . , [u„, t>„], .... We will then have 

«n - Avj; = k (59) 

The sequence of pairs \u u t^], [u 2 , u 2 ], ..., [u„, t>„], ... will 
of course be part of the sequence of numerators and denominators of 
the even convergents of a. If we could assert that k = 1, then we would 
have proved that equation (51) has an infinite number of integral 
solutions. Since we cannot assert this, let us assume that 
k > 1 (in the contrary case when k = 1 everything is proved), 
and go on to the second step of our proof. 

We shall now prove that among the pairs of integers [u 1; Uj], 
[u 2 , t> 2 ], . . . , [u„, u n ], . . . there will be infinitely many pairs 
yielding the same remainders when divided by k. To put it 
another way, we shall prove that there exist two nonnegative 
integers, p and q, both less than k, such that for an infinite 
number of pairs [u u Uj], [u 2 , u 2 ], ..., [ u „ , t>„], ... the equalities 

u„ = a„k + p, v n — b„k + q (60) 

hold, where a„ and b n are the quotients upon division of u„ 
and v n by k, and p and q the remainders. For, if we divide u„ and v„ 
by the integer k, k> 1, then we obtain relations of the form (60), 
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where as always the remainders upon division lie between zero 
and k — 1. Since the only possible remainders upon the division 
of the numbers u„ by k are the numbers 0, 1, 2, .... k — 1, and 
likewise the remainders upon the division of v n by k can only be these 
same numbers 0, 1, 2, . . . , k — 1, then the number of possible pairs 
of remainders upon the division of the numbers u„ and v„ by k 
will be k ■ k = k 2 . This is also obvious because a pair of remainders 
[p„, <?*] corresponds to each pair [u„, u„] and the number of 
different values assumed by each of the numbers p„ and q n separately 
is not greater than k. Consequently, the number of different 
pairs of remainders is not greater than k 2 . Thus to each pair of 
integers [u„, t>„] there corresponds a pair of remainders [p„, <?„] on 
division by k. But the number of different pairs of remainders is 
finite, does not exceed k 2 , while the number of pairs [u„, i>„] is 
infinite. This means that since the number of different pairs in the 
sequence [p,, qj], [p 2 , q 2 \ ..., [p„, <?„], ... is finite, at least 
one pair of remainders is repeated an infinite number of times. 
Denoting this pair of remainders [p, q], we see that there exists an 
infinite set of pairs [u„, r„] for which relations (60) hold. Since 
not all the pairs satisfy (60) for certain definite p and q, whose 
existence we have just proved, we shall renumber all those pairs 
[u„, t>„] which satisfy (60) denoting them by [/?„, S„], And so, the 
infinite sequence of pairs [Rj, S,], [R 2 , S 2 ], . . . , [f?„, S„], ... is 
a subsequence of the sequence [u„, u„] which, in turn, is a subsequence 
of the sequence of numerators and denominators of the even 
convergents of a. The pairs of numbers [/?,, Sj], [/? 2 , S 2 ], ..., 
. . . , [R m S„], . . . satisfy equation (59) and yield the same remainders, 
p and q, on division by k. 

Now that we have established the existence of an infinite set 
of such pairs of positive integers R n and S n , we can go on to the third 
and last step of our proof. Note first of all that the pairs 
[*„, S„l being the numerators and denominators of convergents, 
must be pairs of relatively prime numbers, i. e. pairs of 
numbers which do not possess common divisors. Indeed, if we 


replace k by 2k in relation (24) and set 8 2t = - / ^ L , § 2 *_,= 

TC2k 

P 

= — , then from the equation 

0.2k.- 1 

?2k P 2k- 1 1 


(?2k Q2k-1 Q,2kQ.2k-\ 


multiplying both sides by Q 2k Q 2it ~ ], we get 


p 2kQ2k- 1 — Q2kP 2k- 1 ~ 1 ( 61 ) 
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This relation between four integers, P 2k , Q 2k , P 2k -i 
and Q 2k - u shows that if P 2k and Q 2k have a common divisor greater 
than unity, then its whole left-hand side must be divisible by this 
common divisor. But the right-hand side of equality (61) is unity, 
which cannot be divided by any integer greater than unity. 
Thus it is established that the numbers R n and S n , which can 
only be the numerators and denominators of convergents, are 
relatively prime. From relation (7) it also immediately follows that 

Q 2 <Q 4 <...<Q 2n <... 

From the fact that the numbers R n and S n are relatively prime 
and from the fact that the numbers S lt S 2 , ..., S„, ..., which are 
taken, from the sequence of numbers Q 2n all differing from one 
another, are also all different from one another, it immediately 
follows that in the infinite sequence of fractions 
Kl Rl Rn 

sr ■ s 2 s „’ -• 

there are no numbers equal to one another. Let us write two 
equalities following from the definition of the numbers R n and S„: 

Rl-ASj =(«! -aSiMRj +otS,) = k (62) 

and 

R\ - ASj = ( R 2 - a S 2 )(R 2 + a S 2 ) = k (63) 

where, as above, 'a = ]/~A. 

Also, 

(R i — &S j) (R 2 + aS 2 ) = R\R 2 — /1SiS 2 + cc(RiS 2 Sji? 2 ) (64) 
since a 2 = A. Similarly 

(Ri + ctSi) (R 2 — aS 2 ) = R 2 R 2 — ASjS 2 — tx(i? 1 S 2 S 2 R 2 ) (65) 

When divided by k, R n and S n leave remainders p and q independent 
of n. Consequently, because of relations (60), 

R n = c n k + p, S„ = djk + q (66) 

A series of easy transformations and substitutions leads to 

RiR 2 — v4S^S 2 = Ri ( c 2 k p) — {d 2 k + q) = 

= Ri [(c 2 - Ci)k + Cik + p] - ASj [_(d 2 - djk + d t k + <?] = 

= Ri C(c 2 - c,)k + /?,] - [(d 2 -d,)k + S,] = (67) 

= k[Ri(c 2 — c 2 ) — y4S t (d 2 — dif] + R 2 — AS 2 = 

= k (c 2 — Ci) — (d 2 — di) + lj = kxi 
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where x, is an integer because R 2 — AS\ = fe. Analogously 
RiS 2 — S t R 2 = 

= R, [(d 2 - di)k + d v k + q] - S, [(c 2 - c,)fe + c,fe + p] = 

= R, [(d 2 — d,)fc + S,] — Si l(c 2 — Ci)k + R,] = 

= k [ Ri (d 2 - di) - Si (c 2 - c,)] = fey, 

where y, is again an integer. We can assert that y, is not equal to 
zero. For suppose y, = 0, then 

fey, = RiS 2 — R 2 Si = 0 

whence 

ZL = ^L 
s, s 2 

which is impossible as we established that all the fractions 
RJS n are different. 

Equalities (67) and (65) show that 

(Ri - aS,)(R 2 + aS 2 ) = fex, 4- afey, = fe(x, + ay,) (69) 

and 

(R, + a S,)(R 2 - aS 2 ) = fex, - a fey, = fc(x, - ay,) (70) 

Multiplying (62) and (63) termwise and taking into account 
expressions (69) and (70), we arrive at 

k 2 =(Rl-ASl)(R 2 2 -ASl) = 

= (R, - aS,)(R 2 + aS 2 )(R, + aS,)(R 2 - aS 2 ) = (71) 

= k 2 (x i + ay,)(x, - ay,) = fe 2 (xf - Ayl) 

Cancelling fe 2 out of the result, we finally get 

xj - Ayl = 1 (72) 

But y, 5 ^ 0 which means that x, / 0, otherwise the left-hand side 
would be negative while the right-hand side would be equal to 
unity. Thus, even under the assumption that fe # 1, we have 
determined two non-zero integers, x, and y,, which satisfy equation 
(51). The theory of equations of this type is now complete since 
we know that they do possess a solution for any positive 
integer A and irrational \[A\ and we know how to construct all 

the solutions with the aid of the least solution whose existence 
has been proved. 
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In practice the least solutions should be sought by trial 
and error, choosing the values of x 0 and y 0 . 

We have fully treated the equation 

x 2 — Ay 2 = 1 

when A > 0 and a = ]/~A is irrational. 

If A > 0 and a = ]/~/ 1 is an integer, then this equation may be 
written in the form 

x 2 - a 2 y 2 = (x + a y)(x - ay) = 1 

and since a is an integer and x 0 , y 0 are integers satisfying the 
equation, we must have 

*o + o^o = 1» *o ~ «yo = 1 

or 

x 0 + oty 0 = - 1, x 0 - ay 0 = - 1 

because the product of two integers is unity if and only if each of 
these integers is either +1 or — 1. Each of the systems of two 
equations in two unknowns x 0 and y 0 admits of only trivial 
integral solutions, x 0 = 1, y 0 — 0 and x 0 = — 1, y 0 = 0 respectively. 
Hence, when A is equal to the square of an integer, equation (51) 
has only trivial solutions in integers x 0 = ± 1, y 0 = 0. When A 
is a negative integer, equation (51) has the same trivial integral 
solutions. (When A = — 1, it also has the (symmetrical) trivial 
solutions x 0 = 0, y 0 = ±1.) 

Let us now consider a more general equation 

x 2 — Ay 2 = C (73) 

where A and C are integers, A is positive and a = ]/a is irrational. 
We have already seen that when C = 1 this equation always 
possesses an infinite number of integral solutions (see Theorem III). 
But for arbitrary values of C and A, this equation may not have a 
solution at all. 

EXAMPLE. Show that the equation 

x 2 — 3y 2 = — 1 (74) 

possesses no integral solution. First note that the square of an odd 
number, when divided by 8, always leaves a remainder equal to 1. 
Indeed, since any odd number a may be written as a = 2N + 1, 
where N is an integer, we have 

a 2 = (2 N + l) 2 = 41V 2 + 41V + 1 = 4JV(JV + 1) + 1 = 8M + 1 (75) 
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where M is an integer, as either N or N + 1 must be even. 
Suppose now that [x 0 , y 0 ] is a solution of equation (74). Then 
the two numbers, x 0 and y 0 , cannot have the same parity, that is, one 
must be even, and the other odd. If x 0 and y 0 were both either 
even or odd, then x© — 3 yl would be an even number and 
could not be equal to 1. If x 0 were odd, and y 0 even, 
then Xq divided by 4 would leave a remainder of 1, while 
- 3 yl would be divisible by 4. Hence the remainder upon the division 
of Xq - 3yo by 4 would be unity. This is impossible because the 
right-hand side of equation (74), when divided by 4, leaves a remainder 
of — 1 or 3 = 4—1. Lastly, if x 0 is even and y 0 is odd, then 
Xq is divisible by 4, and, according to (75), - 3 yl may be written 
in the form 

-3y§ = -3(8 Af + 1)= — 24M — 3 = 4(— 6M — 1) + 1 

and so yields a remainder of 1 when divided by 4. Hence 
once again the remainder upon the division of Xq — 3 yl by 4 must 
be 1, which is impossible, as we have seen. Thus, no integers x 0 , 
y 0 can satisfy equation (74). 

We shall not consider the problem of specifying conditions, 
imposed on C and A, under which equation (73) will have 
solutions. It is a difficult one and is solved with the aid of the 
general theory of quadratic irrationalities which belongs to the 
algebraic number theory. We shall only discuss the case when 
equation (73) has non-trivial solutions. As above, we shall call a 
solution [x', y'] non-trivial if x', y' # 0. 

Thus, suppose equation (73) admits of a non-trivial solution 
[x', y'], in other words, suppose that 

x' 2 - Ay' 2 >= C (76) 

Consider for the same A equation 

x 2 Ay 2 = 1 (77) 

For A > 0 and a = J /a irrational this latter equation possesses 
an infinite number of integral solutions [x, y], each of which may be 
determined by 

x = ±x„, y = ±y„ 

where x„ and y„ are determined by formulae (50) of § 4. Since 
[x, y] is a solution of equation (77), 

x 2 - Ay 2 = (x + ay)(x — ay) = 1 

In its turn, equality (76) can be written as 
(x' + a/)(x' - ay') = C 
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( 78 ) 


Multiplying the last two equations term by term we get 
(*' + ay')(x + ay)(x' — a y')(x - ay) — C 
But 

(x' + a y')(x + ay) = x'x + Ay'y + a (x'y + y'x) 
and, similarly, 

(x' — ay')(x - ay) = x'x + Ay'y - a (x'y + y'x) 

Using these two results, we can rewrite (78) in the form 
[x'x 4- Ay'y + a (x'y + y'x)] [x'x + Ay'y - a (x'y + y'x)] = C 
or as 

(x'x + Ay'y) 2 - A (x'y + y'x) 2 = C 

We have thus proved that if [x', y'] is a solution of equation 
(73) then this equation is also satisfied by the pair of numbers 

l>> yl 

x = x'x + Ay'y, y = x'y + y'x (79) 

where [x, y] is an arbitrary solution of equation (77). Therefore 
we have proved that if equation (73) has at least one 
solution then it has an infinite number of solutions. 

We must not assert of course that formulae (79) provide all the 
solutions of equation (73). In the theory of algebraic numbers 
it is proved that all the integral solutions of equation (73) 
may be obtained by taking a certain finite number of solutions, 
depending on A and C, and generating them with the aid of 
formulae (79). When A is negative or equal to the square of an 
integer, equation (73) cannot have more than a finite number of 
solutions. The proof of this proposition is simple and we leave it for 
the reader. The solutions in integers of the most general 
equation of the second degree in two unknowns 

Ax 2 + Bxy + Cy 2 + Dx + Ey + F — 0 (80) 

where A, B, C, D, E, and F are integers, may be reduced, by 
change of variables, to the solution of equations of type (73) with a 
positive or negative A. Hence the behaviour of the solutions, if they 
exist, is the same as that for equations of type (73). 

Summing up what has been proved, we can say that 
equations of the second degree in two unknowns of type (80) may 
either have no integral solutions at all, or have only a finite number of 
solutions, or they may have an infinite number of solutions. In the 
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last case all the solutions can be obtained from a finite number of 
generalized geometric progressions (79). A comparison of the behaviour 
of integral solutions of equations of the second degree in two 
unknowns with the behaviour of the solutions of linear equations 
reveals an extremely important fact. Whereas the solutions of 
linear equations, if they exist, form arithmetic progressions, the 
solutions of equations of the second degree, when they are 
infinite in number, are taken from a finite number of generalized 
geometric progressions. In other words, pairs of integers which 
provide solutions of an equation of the second degree occur much 
less frequently than in the case of linear equations. This is not 
accidental. It turns out that equations in two unknowns of degree 
higher than the second, generally speaking, have only a finite 
number of solutions. Exceptions to this rule are extremely rare. 

§ 6. Equations in Two Unknowns of Degree 
Higher Than the Second 

Equations in two unknowns of degree higher than the second 
almost always, with rare exceptions, have only a finite number of 
solutions in integers x and y. Let us consider first of all equation 

a 0 x n + aiX n ~ x y + a 2 x"~ 2 y 2 + ... + a n y" = c (81) 

where n is an integer greater than two and all the numbers a 0 , 
a u a 2 , . . . , a n and c are integers. 

At the beginning of this century, A. Thue proved that this 
equation possesses only a finite number of solutions in integers x 
and y, with the possible exception of cases when the homogeneous 
left-hand side is a power 

(1) of a homogeneous linear binomial 

(ax + by)" = c 0 

or (2) of a homogeneous quadratic trinomial 
(ax 2 + bxy + cy 2 ) n = c 0 

In both these instances integral solutions can exist only if c 0 is the 
nth power of some integer and, consequently, if equation (81) 
reduces to an equation of the first or of the second degree 
respectively. 

Thue’s method is too complicated for us to describe here. 
We shall confine ourselves to a few notes explaining how the 
finiteness of the number of solutions of equation (81) is 
demonstrated. 
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Divide both sides of equation (81) by y", then 



For the sake of simplicity we shall suppose not only that all the 
roots of the equation 

a 0 z" + a k z"~ l + ... + a n .^ + a„ = 0 (83) 

are different from one another and a 0 a„ + 0, but also that these 
roots cannot be the solutions of any equations of a lower 
degree with integral coefficients. This case is the essential one for 
our discussion. 

In courses of higher algebra it is proved that any algebraic 
equation has at least one root, whence, since any polynomial in z 
is divisible by z — a if a is its root, it follows very simply that 
any polynomial may be represented as 

a 0 z" + cqz" - 1 + ... + a„ = a 0 (z - aj)(z - a 2 )...(z - a„) (84) 

where a,, a 2 , ..., a„ are its n (different) roots. Using this 
expression for a polynomial in the form of a product, we can 
rewrite equation (82) in the form 



Suppose there exist an infinite number of integral solutions 
[x*, y t ] to equation (85). This means that there exist solutions 
with y* arbitrarily large in absolute value. If there existed an 
infinite number of pairs with y k bounded, less in absolute value 
than some definite number, and with x k arbitrarily large, then 
there would be a contradiction, as with such x k the left-hand side 
would be arbitrarily large, while the right-hand side would remain 
bounded. Now suppose y k is a very large number. Then the 
right-hand side of equation (85) will be small and this means 
its left-hand side must also be small. But the left-hand side of 
the equation is a product of n factors containing xjy k and a 0 
which, being an integer, is not less than unity. Consequently, 
the left term can become small only if at least one of the 
factors 



is small in modulus. Clearly, this difference can be small only when 
a m is real, in other words, when the relation a m = a + bi, 
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b 0 does not hold. In the opposite case the modulus of our 
difference cannot be arbitrarily small, since 


x k , . 

a — bi 



yk 



- k -a) +b 2 >\b\ 

Vt / 

Two differences, that is two factors of the left member of equation 
(85), cannot be small in modulus simultaneously because 


Xk_ 

yk 


— a m — 


yk 


- a, 


= |a m -a s 1^0 


(86) 


as the numbers a m are all different. If one of the two factors 
is less in modulus than — | a m — a s |, then, because of relation 
(86), the other one must be greater than -y | a m — a s | . This is 


a consequence of the fact that the absolute value of a sum does 
not exceed the sum of the absolute values. Since the numbers 
<x m are all different, the smallest difference in absolute values, 
| a m - a, | , will be greater than zero (m / s). Denoting it by 2d, 
we will have that if for some sufficiently large y k (which we 
can assume since y k increases indefinitely), 



< d 


then 


Xk_ 

yk 


> d, s = 1, 2, ..., n, s^ m 


(87) 


Then, since the modulus of a product is equal to the product 
of the moduli of its factors, it follows from equation (85) that 


x k 


x * ~ 

x * ~ 

x * ~ 

<*! 


a m- 1 

“m 

a m+ 1 

yk 


yk 

yk 

yk 


yk 




( 88 ) 


If in this equation we replace each of the differences 



, s ^ m by the smaller quantity d and replace |« 0 | by 
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unity, which must be smaller than the integer | a 0 1, then the left-hand 
side of equality (88) will become less than the right-hand side, 
and we get the inequality 


yk 


— a„ 



or 


Xk_ 

yk 


- 




(89) 


where c t does not depend on x„ and y n . There are not more 
than n numbers a m , while the set of pairs [x k , y k ] which satisfy 
inequality (89) for some m is infinite. Therefore, there exists a certain 
m for which, with the corresponding a m , this inequality is valid 
an infinite number of times. In other worlds, if equation (81) has 
an infinite number of integral solutions, then the algebraic equa- 
tion (83) with integral coefficients possesses a root a for which 
inequality 

ot- ' 

q 

holds for arbitrarily large values of q. Here, p and q are integers, 
A is a constant, independent of them, and n is the degree of an 
equation of which a is a root. 

If a were an arbitrary real number, it would have been possible 
to choose it so that there were indeed an infinite number of 
solutions to equation (90) in integers p and q. But in our case 
a is the root of an algebraic equation with integral coefficients. 
Such numbers are called algebraic and they possess special 
properties. The degree of an algebraic number is the degree of the 
algebraic equation of least degree with integral coefficients, which 
is satisfied by this number. 

A. Thue proved that for an algebraic number a of degree n the 
inequality 



<— l — , 3 (91) 

4+ 1 
<r 

can have only a finite number of integral solutions [p, q]. But 
if n > 3 the right-hand side of inequality (90) for a sufficiently 
large q will become less than the right-hand side of inequality 

(91) since n> ~ + 1. Therefore, if inequality (91) can have only a 



47 



finite number of solutions in integers p and q, then inequality 
(90) must certainly have only a finite number of solutions. This 
means that equation (81) can have only a finite number of solutions 
in integers when all the roots of equation (83) cannot be roots of an 
equation with integral coefficients of degree lower than n. It is not 
difficult to ascertain that for n = 2 and given A inequality (90) 
does indeed have an infinite number of integral solutions in p and q. 
The theorem by A. Thue was subsequently strengthened significantly. 
It should be mentioned that the method he used to prove 
his theorem did not allow him to find an upper bound for the 
solutions, in other words, a bound for the possible values 
of | x 1 and | y | for given coefficients a 0 , a u ..., and c. 
This question still remains open. However, the method due to Thue 
enables us to discover an upper bound, though a rather crude one, 
for the number of solutions of equation (83). For certain classes of 
equations of type (83) this bound may be made much more precise. 
For example, the Soviet mathematician B. N. Delone proved 
that, except for the trivial solution [0, 1], equation 

ax 3 + y 3 = 1 

where a is an integer, cannot have more than one integral 
solution [x, y]. He also demonstrated that equation 

ax 3 + bx 2 y + cxy 2 + dy 3 = 1 

where the coefficients a, b, c, d are integers can have no more 
than five solutions. 

Let now P(x, y) denote an arbitrary polynomial in x and y with 
integral coefficients A ks : 

P(x,y) = ^4x‘/ 

We shall say that the polynomial is irreducible if it cannot be 
represented as a product of two other polynomials with integral 
coefficients, each of which is not equal to a number. 

By using a special and extremely complicated method Siegel 
proved that equation 

P(x, y) = 0 

where P(x, y) is an irreducible polynomial in x and y of degree 
higher than the second (i. e. a polynomial which includes terms, 
or a single term, A ks x k y s with k + s > 2), may have an infinite 
number of integral solutions [x, y] only when there exist 
numbers a n , a n _ l5 ..., a 0 , a_„ ..., a_„ and b n , ..., b 0 , 

b- u .... b_„, where n is some integer, such that when expressions 
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where n is some integer, such that when expressions 
x = a„t" + + ... + a 0 + ^j-+... + 

b b 

y — b„t" + b„~ it" - 1 + ... + b 0 -\ + . . . -t — 

are substituted into our equation for x and y an identity 

P(x, y) s 0 

is obtained with respect to t. 

§ 7. Algebraic Equations in Three Unknowns of Degree 

Higher Than the Second. Some Exponential Equations 

Is the number of integral solutions of an equation finite or not? 
Though we can give an answer to this question for equations in 
tw^ unknowns, we can only answer it for very particular 
types of equations in three or more unknowns of degree higher 
than the second. Nevertheless, in these particular cases a more 
difficult problem of actually determining all the integral solutions 
can be solved. Consider for example the so-called Fermat’s 
last theorem. Pierre Fermat, an eminent French mathematician, 
asserted that for any integer n > 3 equation 

x" + f = z" (92) 

has no solutions in positive integers x, y, z. The case xyz = 0 is 
excluded by the requirement that the unknowns be positive. He even 
claimed to have a proof of this proposition (evidently, using 
the method of infinite descent, see below), but it has never been 
found. When the German mathematician E. Kummer subsequently 
attempted to prove Fermat’s theorem, he for some time thought he 
has succeeded. However he discovered that one proposition, correct 
for usual integers, does not hold for the more complicated number 
formations which naturally arise in research connected with 
the problem. This was that the factorization of what are called 
algebraic integers, in other words, roots of algebraic equations with 
integral rational coefficients and with a unit coefficient of the 
leading term, into undecomposable prime integral factors of the same 
algebraic nature is not unique. The factorization of the ordinary 
integers is of course unique. For example, 6 = 2-3, with no 
other factorization being possible in the domain of ordinary integers. 

Consider now the set of all algebraic integers of the type 
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m + n J/—5 where m and n are ordinary integers, and note that 
both the sum and the product of two such numbers are again 
numbers of the same set. A set of numbers which contains any 
sums and products of the numbers in it is called a ring. 
By definition the ring under discussion contains the numbers 2, 3, 

i + ]/-~5 and 1 — \/—5. It is easy to ascertain that each of these 
numbers is a prime; none of them can be represented as a product 
of two integers of the ring neither of which are equal to unity. However, 

6 = 2-3 = (1 + 1/^5) (1 - ]/^S) 

so that in our ring the number 6 does not factorize uniquely into 
prime factors. 

Non-uniqueness in the factorization into prime factors also 
occurs in other, more complicated, rings of algebraic integers. 
Having discovered this, Kummer realized that his proof of Fermat’s 
theorem in the general case was false. In order to overcome the 
difficulties connected with the non-uniqueness of the factorization 
into factors Kummer constructed the theory of ideals, which is 
extremely important in modern algebra and number theory. 
Even with the aid of his new theory, Kummer was unable to prove 
Fermat’s theorem in the general case, and proved it only for 
those n which are divisible by at least one of what are known as regular 
prime numbers. We shall not go into what is meant by the 
concept of regular prime number. It is not even known at the present 
time whether there is a finite number of regular prime numbers 
or infinitely many of them. 

At the present time, Fermat’s last theorem has been proved for 
many values of n and, in particular, for any n divisible by a prime 
number less than 100. Fermat’s last theorem turned out to be 
extremely important for the development of mathematics in general 
because the attempts to prove it led to the discovery of the 
theory of ideals. It should be mentioned that this theory was 
independently constructed in quite another way and for a different 
reason by E. I. Zolotarev, an eminent Russian mathematician who 
regrettably died in the prime of his creative life. These days, a proof 
of Fermat’s last theorem, especially a proof based on the concepts 
of the theory of divisibility of numbers, would have only curiosity 
value. If, however, a proof were attained by a new and fruitful 
method, then its importance, or, rather, the importance of 
the method itself, might be quite great. Even now amateurs 
continue to attack Fermat’s theorem by elementary methods. 
All such attempts are doomed to failure. Elementary arguments 
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proceeding from the theory of divisibility of numbers were 
used by Kummer and were further developed by some of the most 
eminent mathematicians* yet nothing important was obtained. 

We shall now prove Fermat’s theorem for the case n = 4, 
since the method of infinite descent on which the proof is based is very 
interesting. 

THEOREM IV. Fermat's equation 

x* + y* = z 4 (93) 

has no solutions in integers x, y, z, xyz # 0. 

Proof. We shall prove an even stronger proposition, namely, that 
equation 

x 4 + / = z 2 (94) 

has no solution in integers x, y, z, xyz / 0. From this theorem it 
immediately follows that equation (93) has no solution. If equation 
(94) has a solution in non-zero integers x, y, z, then we may 
assume that these numbers are pairwise relatively prime. For if there 
is a solution in which x and y have a greatest common 
divisor d > 1, then 

x=dx u y = dy t 


where (x 1; y x ) = 1. Dividing both sides of equation (94) by d A , we have 

4 + yi = (jrj = (95) 

But x, and >>, are integers, therefore z, = z/d 2 is also an integer. 
Now, if Zj and had a common divisor k > 1, then, because of 
expression (95), x 2 would have to be divisible by k, which means x, 
and k could not have been relatively prime. Hence we have proved 
that if there exists a solution to equation (94) in non-zero integers, 
then there also exists a solution in non-zero and pairwise 
relatively prime integers. Therefore it is sufficient for us to 
prove that equation (94) does not have solutions in non-zero 
pairwise relatively prime integers. In the following proof, when we say 
that equation (94) has a solution, we mean that it has a solution 
in positive pairwise relatively prime integers. 

In § 3 we proved that all the positive integral and pairwise 
relatively prime solutions of equation (12) 

x 2 + y z = z 2 (96) 


are determined by formula (18) and have the form 

u 2 — v 2 u 2 + v 2 

x = uv, y = r , z = 


(97) 
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where u and v are any pair of odd and relatively prime positive 
numbers. 

Let us give another form for formulae (97) determining all the 
solutions of equation (96). Since u and v are odd numbers, 
then setting 


u + v 
2 


= a, 



(98) 


we determine u and v by 

u = a + b, v = a — b 


(99) 


where a and b are integers with different parity (one is even 
and the other odd). Equalities (98) and (99) show that to any pair of 
odd and relatively prime numbers u and v there corresponds a pair of 
relatively prime numbers a and b of different parity and that 
to any pair of relatively prime numbers a and b of different 
parity there corresponds a pair of relatively prime odd numbers u 
and v. Therefore substituting a and b for u and v respectively 
in formulae (97) we find that all the triplets of positive and pairwise 
relatively prime integers x, y, z (x odd), which are solutions of 
equation (96), are determined by formulae 

x — a 2 -b 2 , y = 2ab, z — a 2 + b 2 (100) 

where a and b are any two relatively prime numbers of different 
parity, on the condition that x > 0. These formulae show that the 
two numbers, x and y, are of different parity. Now, if 
[x 0 , y 0 , 2 0 ] is a solution of equation (94), then 

[4] 2 + [y§] 2 = zl 

so that the triplet [ xl , yl, z 0 ] satisfies equation (96). But then 
there must exist two relatively prime numbers a and b, a > b, of 
different parity, such that 

Xo — a 2 — b 2 , yo = 2ab, z 0 = a 2 + b 2 (101) 

We have assumed here for the sake of definiteness that x 0 is 

odd and y 0 even. In the contrary case nothing is changed since x 0 
could be changed for y 0 and vice versa. We already known from 

equality (75) of § 5 that the square of an odd number 

divided by 4 leaves a remainder of 1. Therefore from equality 

x 2 0 =a 2 -b 2 (102) 

it follows that a is odd and b even. If otherwise, the left-hand 
side of equality (102) divided by 4 would leave a remainder of 1 
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while the right-hand side would leave a remainder of — 1, as we 
assumed a to be even and b odd. Since a is odd and (a, b) — 1, 
we have (a, 2b) = 1. But then froln equality 

yl = 2. ba 

it follows that 


a = t 2 , 2 b = s 2 


(103) 


where t and s are some integers. But it follows from relation 
(102) that [x 0 , b, a] is a solution of equation (96) and therefore 

x 0 = m 2 — n 2 , b = 2 mn, a = m 2 + n 2 


where m and n are some relatively prime numbers of different parity. 
From (103) we have 



whence, as m and n are relatively prime, it follows that 

m = p 2 , n = q 2 (104) 


where p and q are non-zero integers. Since a = t 2 and a = m 2 + n 2 , 
it follows that 


q* + p* = t 2 


(105) 


But 


z 0 = a 2 + b 2 > a 2 


Therefore 

0 < t = ]/a < fa < z 0 (z 0 < 1) (106) 

Setting q = x u p — yy and t = z t we see that if there exists a 
solution [x 0 , y 0 , z 0 ], then there must exist another solution 
[x„ y t , Zj] for which 0 < z x < z 0 . This process of obtaining 
solutions of equation (94) may be continued indefinitely, and we 
obtain a sequence of solutions 

[x 0 , y 0 , Zq], [xj, y u z t ], .... [x„, y m z„], ... 

in which the positive integers z 0 , z u z 2 , ..., z„, ... decrease 
monotonically; in other words, inequalities 

Z 0 > Zy> Z 2 > ■■■> z n > ... 

hold for them. But positive integers cannot form an infinite and 
monotonically decreasing sequence as there cannot be more than z 0 
terms in it. We have thus come to a contradiction by assuming 
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that equation (94) has at least one solution in integers x, y, z, xyz t- 0. 
This serves as a proof that equation (94) does not have a solution. 
Accordingly equation (93) has no solutions in positive integers 
[x, y, z] either, since, if otherwise, if [x, y, z] were a solution 
of equation (93), then [x, y, z z ] would be a solution to (94). 

The method of proof we have employed, consisting in using 
one solution to construct an innumerable sequence of solutions with 
indefinitely decreasing positive z, is called the method of infinite 
descent. 

As was remarked above, Fermat’s last theorem in the general 
case does not yet yield to this method because of the non- 
uniqueness of the factorization of the integers of an algebraic 
ring into prime factors from the same ring. 

Note that we have demonstrated the non-existence of integral 
solutions not only for equation (94), but also for equation 

x 4 " + y 4 " = z 2 " 

It is interesting to note that equation 

x 4 + y 2 = z 2 

possesses an infinite number of positive integral solutions. For 
example, one solution is [2, 3, 5]. It is left to the reader to find 
general expressions for all such solutions of this equation in 
integers x, y, z. 

We shall now consider another example which illustrates the 
method of infinite descent, but the argument will be somewhat 
different. 

EXAMPLE. Prove that equation 

x 4 + 2 y 4 = z 2 (107) 

has no solutions in non-zero integers x, y, z. Let us suppose that 
a positive integral solution [x 0 , y 0 , Zo] does exist. These 
numbers may be immediately assumed to be relatively prime, since 
if they had a greatest common divisor d > 1 , then the numbers 

-T-, —- 7 —, ~~ would also be solutions of equation (107). 
d d d 

Moreover, the existence of a common divisor for any two numbers 
would mean that all three of them had a common divisor. 
Let us also assume that z 0 is the least of all the possible 
values of z in the positive integral solutions of equation (107). 
Now since [x 0 , y 0 , z o] satisfies equation (107), [xq, yl, z 0 ] will be a 
solution of equation 

x 2 + 2 y 2 = z 2 (108) 
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Using formulae (19') from § 3 which provide ail the positive 
integral solutions of equation (108), we see that there exist positive 
integers a and b with (a, b ) .= 1 and a odd, such that 

Xq = ± ( a 2 - 2b 2 ), yl = lab, z 0 = a 2 + 2 b 2 (109) 


From yl = lab it follows that b must be even since y 0 is even, 
yl is divisible by 4 and a is odd. Now as b/2 and a are relatively 
prime, equality 


t-J-4 


directly yields 



where m and n are positive integers and (m. In) = 1. But from 
equalities (109) it follows that 

xl= ±(a 2 -2b 2 )= ±^a 2 j (110) 

where x 0 and a are odd. We have seen that the square of an odd 
number divided by 4 leaves a remainder of 1. Therefore, the 
left-hand side of equality (1 10), upon division by 4, gives a remainder 

of 1 while a 2 — 8 ; when divided by 4, also leaves a remainder 

of 1. This means that the bracket on the right-hand side of (110) 
may be taken only with the plus sign. Now (110) may be written 
in the form 

Xq = m 4 — 8 n 4 


or in the form 

xl + 2(2n 2 ) 2 = {m 2 ) 2 (111) 

where x 0 , n and m are positive and relatively prime integers. 
Thus, the numbers x 0 . In 2 and m 2 constitute a solution of equation 
(108) and are relatively prime. Again, according to formula (19') 
of § 3, there may be found integers p and q, with p odd and 
(p, q) — 1, such that 

In 2 - 2pq, m 2 = p 2 + lq 2 , x 0 — ± (p 2 — lq 2 ) (112) 
But since (p, q) = 1 and n 2 = pq, we have 

p = s 2 , q = r 2 

where s and r are relatively prime integers. Finally, from here 
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there follows the relation 


s 4 + 2 r 4 = m 2 (113) 

which shows that the triplet s, r, m is a solution of equation 
(107). But from the results 

z 0 = a 2 4- 2 b\ a = m 2 

obtained above it follows that z 0 > m. Thus, proceeding from a 
solution [x 0 , y 0 , z 0 ], we haye found another solution [s, r, m], 
in which 0 < m < z 0 . This contradicts the assumption we made 
that z 0 was the least possible value. Thus we have arrived at a 
contradiction by assuming the existence of a solution to equation 
(107), and we have proved that this equation is unsolvable in non-zero 
integers. 

We leave it to the reader to show that equations 
x 4- + 4y 4 = z 2 , x 4 — y* = z 2 
x 4 — y* = 2 z 2 , x 4 — 4y 4 = z 2 

have no positive integral solutions. 

We shall conclude with a few remarks about exponential 
equations. The equation 

a x + b 3 ' = c z . (114) 

where a, b and c are integers, not equal to any power of 2 or to zero 
can have no more than a finite number of solutions in integers 
x, y, z. The same proposition with a weak condition being added is 
valid for arbitrary algebraic numbers a, b and c. Moreover, equation 

A a? ...«?■ + B \ . . . pj- + Cy\ l = 0 (115) 

where A, B and C, ABC # 0, are integers, a lt a„, p 1; . . . 
. . . , p„, and Yj, . . . , y„ are integers and the numbers 

a = <*!... a B , P=P t ... P m , y = Yi ••• Y P 

are relatively prime, can only have a finite number of integral 
solutions [xx, ..., x„ y u .... y m , z lf .... z p ]. A generalization 
of this proposition with A, B and C and a ( , p k and y s being 
algebraic numbers is also possible. Equations of the type (115) and 
their generalizations are of great interest because, as is shown in the 
theory of algebraic numbers, to each algebraic equation of type 
(81), there corresponds a certain exponential equation of type (115) 
and to each solution of equation (81) there corresponds a solution 
of equation (115) in integers. This correspondence extends to equa- 
tions of a more general type than (81) and (115). 
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The book is devoted to one of the most interesting branches of 
number theory, the solution of equations in integers. The solution 
in integers of algebraic equations in more than one unknown with 
integral coefficients is a most difficult problem in the theory of 
numbers. The theoretical importance of equations with integral 
coefficients is quite great as they are closely connected with many 
problems of number theory. Moreover, these equations are 
sometimes encountered in in physics and so they are also 
important in practice. The elements of the theory of equations with 
integral coefficients as presented in this book are suitable for 
broadening the mathematical outlook of high-school students and 
students of pedagogical institutes. Some of the main results in the 
theory of the solution of equations in integers have been given and 
proofs of the theorems involved are supplied when they are 
sufficiently simple. 
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